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o '• 

Can local fluctuations of a "Quintessence" scalar field play a dynamical role in the gravitational 
' clustering and cosmic structure formation process? We address this question in the general frame- 

work of scalar-tensor theories of gravity. Non-linear energy density perturbations, both in the scalar 
field and matter component, and linear metric perturbations are accounted for in the perturbed Ein- 
stein's equations. We derive the Newtonian limit of the relevant equations for clustering in scalar- 
tensor cosmologies. We then specialize to non-linear perturbations of the "Extended Quintessence" 
j — ■ model of Dark Energy; in such a model, a non-minimally coupled scalar field is thought to be re- 

' sponsible for driving the present accelerated phase of the Universe expansion. The interplay between 

Dark Energy and Dark Matter is displayed in the equations governing the growth of structure in 
£SJ , the Universe. 

> ' 
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V/-) , I. INTRODUCTION 

m ■ 

■ Our traditional picture of the Universe has been definitely upset when, in 1998, astronomers found that distant type 
la Supernovae were dimmer than expected in a decelerating Universe Q, Q. This early evidence has been confirmed 
by the subsequent studies, which combined the most recent Cosmic Microwave Background (CMB) data from WMAP 
i-£h ■ (see and references therein) and Large Scale Structure (LSS) data 0, 0,01' together with measurements of the 
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P" 1 ' Hubble constant |8j; the fact that the Universe turns out to be geometrically close to flat, together with the estimates 
of its matter content, has called for deep changes of the old "standard" scenario of a matter dominated Universe. Ac- 
cording to these observations, almost 70% of the total energy of the Universe resides in a "Dark Energy" component, 
which is plausibly acting as a repulsive force driving the cosmic acceleration. Although the Cosmological Constant 
has been historically proposed as a candidate for the cosmic acceleration, the theoretical difficulties in justifying its 
, exceedingly small value motivated the search for alternative theories. 

There are currently two different general approaches to Dark Energy modeling; one class of models is based on the 
introduction of a new cosmological component having negative equation of state, often described through the dynam- 



ics of a self-interacting, minimally coupled scalar field 9]- 19] . The other class of models proposes modifications to 
the gravity itself, introducing a non-minimally coupled scalar field (scalar-tensor theories) or changing the function of 
the Ricci scalar appearing in the Gravitational sector of the Lagrangian of the theory j^-liil- I n both pictures, an 
open important issue is whether modifications of gravity, or small-scale perturbations of a Quintessence scalar field, 
could lead to significant effects on the formation of structures, such as galaxies and clusters (see [2^. |2!|,|34|-0])- 
A powerful tool for the investigation of the Dark Energy main properties and parameters has been the linear per- 
turbation theory, which allowed to make accurate predictions on many cosmological observables, such as the CMB 
spectrum of anisotropies (see, e.g., |43j]-|47j): in this paper, we deal with the behavior of perturbations in the non- 
linear regime. We focus on scalar-tensor theories, where a scalar field non-minimally coupled to the Ricci scalar is 
proposed as a candidate for the Dark Energy component; in principle, the "Extended Quintessence" field l21l allows 
for non- vanishin g sm all-scale perturbations, which have been analyzed in the linear regime in |2l|, [2? | . |29 | . [ill ]. 

In particular, in Tisl and E|l it was shown that, while perturbations in a minimally-coupled scalar field behave as 
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radiation on sub- horizon scales 49], so that the field rapidly becomes a smooth component on such scales, perturba- 
tions in a non-minimally coupled scalar field can be dragged by perturbations in the matter component, thanks to 
the coupling of the scalar field to the Ricci scalar, and become non-linear. In practice, the mechanism which would 
damp out the field perturbations is here counterbalanced by the presence of a source term, which is directly related 



2 



to the gravitational potentials and, ultimately, to the matter perturbations: we expect this dragging to be able to 
affect gravitational clustering down to galactic scales. 

We extend here the analysis performed in [2^ for the linear perturbations in scalar-tensor theories, obtaining the 
non-linear Newtonian limit of Einstein's equations for this class of models. These will be the equations to be eventu- 
ally inserted in a modified N-body code, in order to simulate the behavior of collapsing matter under this modified 
property o f gr avity. Some numerical simulations, involving Dark Energy in several context, have been performed in 
Refs. 0-0. 

Finally, we focus on the Extended Quintessence model, specializing the Poisson equation for the gravitational poten- 
tial in order to evaluate modifications with respect to the "standard" theory. 

The plan of the paper is as follows: in section [HJ we write the perturbed Einstein's equations and the equation for 
scalar field perturbations in non-linear regime, assuming that the metric perturbations are linear, while those of matter 
and scalar field are not. In sect. IHII we discuss the approximations giving the Newtonian limit of these equations, and 
in sect. II VI we specify the general equations for scalar-tensor theories to the "Extended Quintessence" case. Finally, 
in sect. [V]we draw our conclusions. 



II. PERTURBED EINSTEIN EQUATIONS IN EXTENDED QUINTESSENCE 



Since we are interested in the problem of how Dark Energy perturbations could affect the structure formation 
process, excluding the limit of very strong gravitational fields (typical, for example, of black holes), we will consider 
linear perturbations of the metric tensor, while keeping non-linearity in the matter and scalar field perturbations |40|. 
This also means that we are restricting our study to non-relativistic Dark Matter, whose typical velocity is much 
smaller than the speed of light even in the presence of highly non-linear matter overdensities. 

Metric perturbations can be decomposed into scalar modes, as well as by vector and tensor ones; by taking the spatial 
covariant divergence of the (0-i) Einstein equations and the trace of the (i-j) equations, we will single out scalar 
modes, so that, without loss of generality, the line-element can be written, in conformal Newtonian-gauge formalism, 
in terms of only two gravitational potentials <I> and 4" as 

ds 2 = a 2 (r/)[-(l + 2<P)drj 2 + (1 - 2^)8\dx l dx 1 ] = g^dx^dx" , (1) 

where 77 denotes conformal time; note that we adopt the signature ( H — | — h). The Lagrangian for scalar-tensor theories 
reads 

C = YJ { ^ R) ~ \»^^ C ~ + Matter (2) 

where k = 8ttG* and G* is the bare gravitational constant. 

We will consider a Universe filled by dark matter and a non-minimally coupled scalar field; the matter component will 
be described in terms a discrete set of particles with constant mass m and coordinates x (tj) (a = 1, 2, ...). Following 
[5^ . if we denote the matter four- velocity by = dx^/dx° (x° = 77), the matter stress-energy tensor reads 

Ty ~ a- 2 Pm u»u v (3) 

where we neglected corrections C|u| 2 , and we defined 

p m = ma- 3 ^V 3) (x-x Q ) . (4) 



The scalar field stress-energy tensor reads 



Tf = w 



Oil/ 1 1 OLV 1 1.(7 

9 <P,u<p,0 - tj5 9upP,<T<P' 



-9 av 9»pV+ l -g a »g„ U - RF ^j +g av F tU . fi -g av g vfi F>' r . tr + (± - F) G% (5) 
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where, R is the Ricci scalar and F = ^§r- Restricting ourselves to the simplest case of f(<j>, R) linear in R, (as is the 
case of the "Extended Quintessence" model), it turns out that FR = ^^' R ^ . We also assume that uj — const — 1, so 
that the tensor JSJ becomes 



T a p = 4>,u<t>,p - ~9at34>,cr4> ,a ' - Vg a fj + F Q;/ 3 - gapF^.a + I F \ G a fj 



(6) 

Both the scalar field and the function F((f>(x,rj)) can be written as the superposition of a background component, 
only dependent on time, plus a generally non-linear perturbation, which depends both on the spatial coordinates and 
on time: 



c,77) = 0b(77)+<ty(x,J7) F(0(x,7 ? )) = ^ b (7 ? )+,5F(x,77) , 



(7) 



The evolution of the background quantities is fully determined by the Friedmann equations and the unperturbed 
Klein-Gordon equation, which read, respectively, 



TL 



H 



1 



a 2 (p„ 



>l + a 2 (V)-3HF h 



]=r (a 2 (Pm> + 24>l ~ 2a 2 (V) + 3F b ) 



6F b 



(8) 



(9) 



where (V) = V(4>h{r])) and (p m ) is the background value of the matter energy-density. Finally, the background 
Klein-Gordon equation in scalar-tensor theories reads 

^ + 2H^-^a 2 R h F^ + a 2 V'{^) = Q , (10) 
where a prime denotes differentiation with respect to <j). 

In our approach to the perturbed Einstein's equations, the Ricci and Einstein's tensors will be linearly perturbed, 
while the perturbations in the components of the stress-energy tensors © and © are obtained by subtracting their 
background (mean) values from the fully non-linear ones: 



= -(2?)) . 

As for the (0-0) component ("energy constraint"), we get the following perturbation equation: 



(11) 



V 2 $ + 3H($ + I-) + 3* = 8ttG* 



\ a /a 



+ (^) 2 + 2</> b 



a 2 (V- (V)) - 2a 2 $V + 



+ ^SF ~ ^F($ + <j) + ^V<5F • (Vf - 3V$) - ^\7 2 SF{1 + 2$ + 2*) 



3H5F 
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- F (3H* + 3* + 37i$ + V^$) 



(12) 



where <$gj = (5^(x-x a ). 

Similarly, the equation for the divergence of the (0-«) or "momentum constraint" equation (herefrom latin indices will 
be used to label spatial coordinates) reads 



H V 2 $ + V 2 * = 4ttG 



rnaoT 1 (^]w( a) ) + 20V* • S/S(f> + (1 + 24')</>V 2 ( 5< ? (> + (1 + 2*)V(50 ■ Vc50- 



WSF ■ ( 2ftW - V* ) + W 2 SF ( -H - 2H^ + * ) + WSF ■ (2W - V$) + V 2 <5F(1 + 2*) - FV 2 <S> 

K 



2V5F ■ (7A7$ + V*) + 2 



F I (HV 2 $ + V 2 4- 



(13) 
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Finally, the trace of the {i-j) component is: 



27i$ - H(<$> + 5*) - 4W 2 $ - f + ^V 2 * - ^V 2 $ 



= 8ttG* 




ma o) ; — a 

(a) 



(Pm) J + ~(1 + 24<) W| 2 + a 2 (^ - (V))+ 



F 



-(1 - 2$)<5F - + F(4W$ + + + V<SF 



-2$(W + 2ft 2 ) - W$ - -V 2 $ - - * + -V 2 * 



1 



-V* + -V$) + -V 2 SF(l + 24») 
6 2 6 

-SF{H + 2H 2 



(14) 



The perturbed Klein-Gordon equation gives 

5<j> +2HS(j) - ($ + 3*)(0 b + <50) - V($ - *) • - (1 + 2$ + 2^>)V 2 8<j) + a 2 [^'(0(x, 77)) - V'((f> h )} + 

a 2 .... 

2 a 2 W(0(x, 77)) - —SF'R b + (i^ + <5F')(3W$ + V 2 $ + 3* + 9ft* - 2V 2 *) = 



(15) 



where i?b — 6a 2 (ft + H 2 ) is the Ricci scalar in the unperturbed metric. Finally, the equation describing the motion 
of Dark Matter particles, in the weak-field limit, is given by 



ft - $ - 2* 



(16) 



, x l a 's, which is redundant in the number 



The system Q12 fl -H16 fl is a closed system of equations in the variables $,'3/, 
of equations, as a consequence of the underlying gauge-invariance of Einstein's theory; as we will see, it will be more 
convenient to disregard the momentum constraint, working with the remaining equations. 

All these equations apply to cosmological scales, under the assumption of linear perturbations of the gravitational 
potentials. In the next section, we will determine the system adequate to determine the matter particle motion on the 
scales relevant for structure formation, extending the standard Newtonian approximation to the case of scalar-tensor 
theories. 



III. THE EXTENDED NEWTONIAN APPROXIMATION 

There are many considerations that can help to simplify the equations 112l) - (|16fl . First of all, note that, in the 
weak-field limit, |$| and |^| are of the order of \u z \ 2 -C 1, \u l \ being the typical velocity of non-relativistic matter; 
furthermore, since the characteristic evolution time of $ and \& is Td yn ~ t = ft , we can neglect $, \& with respect 
to ft in the energy constraint equation, in the perturbed Klein-Gordon and in the equation of motion. 
Focusing on the perturbation growth on scales well below the Hubble radius, we can apply further approximations to 
the system. On those scales, $, ft<£> ~ ft 2 $ <C V 2 $ and i', ft* ~ ft 2 * <C V 2 * and we can neglect such terms in the 
energy constraint, Einstein spatial trace and Klein-Gordon equations. We also assume that ft ~ ft 2 , so the last term 
on the right-hand side of the trace equation is negligible, on sub-horizon scales, with respect to the term containing 
V 2 5F. Equations l|12|) - (|16fl reduce therefore to the following system: 



V 2 $ = 8ttG* 



^SF - ~F($ + *) - \^ 2 5F + 3HSF 



+ (^) 2 + 2<h,6j> - a 2 (V - (V)) + 



F ] V z * 

K 



(17) 
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H V 2 $ + V 2 * = 4ttG* 



53 m^- 1 (<ygju| o) ) . + 0V 2 ^ + V<ty • VSfa 



HV 2 8F + W 2 8F - FV 2 <S> + 2 



HV 2 $ + V 2 * 



(18) 
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_ V 2 * - -V 2 $ = 8ttG* 



ma- 1 8$ - a 2 ( Pm ) + i|V50| 2 + a 2 (V - (V))+ 



$F b - IsP - 2TL8F + F(4H$ + -*) + ^V 2 ^ 

2 2 2 6 



1 \ V 2 



(19) 



2W«S<£ - ($ + 3*)0 b - V 2 ^ + a 2 [y'(0(x, ry)) - V'(fa)\ 
a 2 

— J J F 1 ' J R b + (i ;l b + ^')(V 2 $-2V 2 *) = ; 



(20) 



Hi!, 



(21) 



In order to further simplify the equations, we have to estimate the characteristic evolution time for scalar field per- 
turbations. We will follow the approach of 01 > considering the scalar field perturbation as the sum of a "relativistic" 
perturbation 8 fa and a "non-relativistic" one, <5</>nr, 



<50(x, 77) = 8fa(x, 77) + <50 N r(x, rf) 



(22) 



where 8fa^ is defined as the solution of the perturbed Klein-Gordon equation in the c — > 00 limit. By definition, the 
time derivative of <5</>nr is negligible with respect to that of 8 fa, which, on the contrary, has a wave-like behavior. 
The characteristic (conformal) time for the variation of 50n R , tnr, is much larger than the characteristic time for the 
variation of 8 fa, t r : on scales much smaller than the horizon, the time-variation of the non-relativistic component 
of scalar field perturbations is negligible with respect to the spatial gradient, which is equivalent to the limit c — ► 00. 
The Klein- Gordon equation for 50nr reduces to 

- V 2 8fa R + a 2 8V'(6<fa R (x, r))) + (F b + <5i^ R )V 2 ($ - 2*) - y<5F^ R i? b = (23) 

where SV'(8fa R (x, V)) = V'(fa(v) + ^nr(x, r))) - V'(fa(r))), and 8F NR = F(fa + 8 fan) - F(fa). 
Inserting the definition 1221) into eq. H20|) . an d using (|23() . we obtain the equation for 8 fa: 

8 fa + 2H8fa - V 2 8fa, + a 2 [V'(fa + 8 fa, + Stfan) - V'(fa + ^nr)] - y [F'(fa + Stfaa + 5 fa) - 
- F' (fa + 8fa R )]R h + [F^fa + 8fa K + 8fa)-F' (fa + 8 fa K )]V 2 ($-2^ = . (24) 

In the last equation we neglected all terms of order H 2 : since Sfa R J Sfa R ~ t nrj we have 

4r/%r~^<H 2 ; 2H8fa R /8fa R ~ < U 2 ; (<i> + 3V)fa ~ H 2 . (25) 

We assume that, on the scales characteristic for structure formation, 8fa, is a small perturbation (this assumption 
will be verified a posteriori); in such a case, V'(fa + <^>r + <^>Nr) — V(fa + <5</>nr) ~ V"(fa + 8fa R )8fa , and 
F'(fa + 8fa + Sfa R ) — F'(fa + Sfa R ) ~ F"(fa + 8fa R )8fa. Substituting in eq. pijl. there will be two terms 
containing F" (fa + 8fa R )8fa; they will produce, respectively, a term of order 8fak 2 (& — 2^) and a term of order 
8 faH 2 , which are both negligible with respect to V 2 8fa. Therefore, eq. I|24|l becomes 



b R + 2H8fa, - V 



+ a z V 



»NR, 







(26) 



G 



This equation describes a (quasi-massless) plane-wave, whose amplitude is damped by the cosmic expansion, so that 
in a timescale r ~ 7^ _1 the relativistic perturbation 5</>r vanishes on scales smaller that the horizon, — > 0. The 
mass associated to this wave is a 2 V"(4>b + £0nr) ~ H 2 , negligible on sub- horizon scales: as the perturbation £</>r 
enters the horizon, it will behave as radiation. We have thus verified a posteriori the validity of our assumption 
on the smallness of <5</>r, which will be neglected hereafter. From now on, 5(f> has to be understood as <5</>nr, the 
non-relativistic component of the scalar field perturbation, and we will omit the subscript nr.; this will be, in general, 
a non-linear perturbation. 

In the Newtonian limit we can neglect time derivatives in the trace equation and in the energy constraint. Indeed, 
the scales which we are considering are well below the horizon, so the Laplacian terms will dominate over terms of 
order Ti 2 . Note that </> 2 also is of order H 2 , as a consequence of the Friedmann equations and of the smallness of the 
coupling at low redshifts; we will neglect <fi in this limit. 
The energy constraint reads: 



FV 2 $ 



\ ($>"W( ( a) - a 2 ( Pm ) J - a 2 (V- (V)) - \V 2 5F 



(27) 



and the trace equation reads 

^(4V 2 *- V 2 $) 



\ ^ma- 1 6$-a 2 ( Pm )^+±\VS0\ 2 +a 2 (V-(V)) + 



-V 2 5F 
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(28) 



In the equations above, S<f> = 50nr- Note that, as anticipated in the previous section, the system i|17 |l -(|19 |l provides a 
redundant set of equations for the gravitational potentials $ and "J: the (0-0) equation and the trace equation fully 
specify their evolution, so we can get rid of the spatial covariant derivative of the (0-i) Einstein's equations. Using 
eq. lf2*7jl and J2EJ), together with eq. J53J), we solve for 6<fi: 



V 2 * 



1 

2F 



a 2 6p m + a 2 SV + V 2 SF + -|V<5</>| 5 



(29) 



V 2 $ = — [a 2 S Pm - 2a 2 SV - V 2 <5F] ; (30) 
At 

where SV = V(<pb + Sip) — V(</>t>) and a 2 Sp m = i\ ja a^ 1 m a 5^J ) — a 2 (p m ). The perturbed Klein-Gordon equation 
becomes 

F! -\-/)F' n 2 
V 2 50 = a 2 SV - — [a 2 S Pm + 4a 2 6V + 3V 2 5F + |V^| 2 1 - — 5F'R h (31) 

2.T 2 

Once the function F is specified together with the scalar field potential V, the system of equations H21|) . 130fl . 
(|31|l fully determines the evolution of perturbations in matter, scalar field and gravitational potentials; it can be 
numerically integrated, given the appropriate initial and boundary conditions. It is interesting to note that, for the 
"ordinary Quintessence" models, i.e. minimally coupled scalar field, the non-relativistic perturbations of the field 
vanish on sub- horizon scales as a consequence of the vanishing anisotropic stress (see Appendix iBl and : however, 
this is generally not true for the extended models we are dealing with: the scalar field coupling to the Ricci scalar 
affects the final Poisson equation for the gravitational potentials, and can give rise to substantial modifications of the 
standard structure formation picture. 

In the next section, we will consider the special case of non- minimal coupling analyzed in 21], adopting the current, 
local constraints on the value of the coupling function F. 
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IV. EXTENDED QUINTESSENCE 

Now, let us focus on the "Extended Quintessence" model of ; in this model, 

F = k- 1 + C<f> 2 (32) 

where k is a universal constant, whose unknown value is proportional to the "bare" gravitational constant G*; in scalar- 
tensor theories, G* can in principle deviate from the Newtonian constant measured in Cavendish- type experiments [28| . 
Indeed, the actual Newtonian force between two close test masses measured at the present time in such experiments 
is proportional to 



Ge//|o = G = ¥r( UJBD ^ I0 ) ~ ~~rr = Gn ■ .:«! 
nF \lujbd,o + 3/2J kF 

In the equation above, ujjbd,o is the present value of the Jordan-Brans-Dicke (JBD) parameter, which reduces to 

1 F' 2 

- % (34) 



when the kinetic factor us(4>) in the Lagrangian J2J) is chosen to be a constant, as in our case. 

The term between brackets in eq. (I33|l is due to the exchange of a scalar particle between the two test masses: since 
the JBD parameter is bound from solar-system experiments to be larger than - 3000 (0, [6o|). the present values 
of G e ff and Gjv almost coincide. F is the coupling function entering in the Lagrangian of the theory, and, as we 
have seen, it enters in the Friedmann equations and in the Poisson equations for the local gravitational potentials; 
the value of this function at present time can be determined only locally, i.e. on the length scales where we are able 
to perform gravitational experiments to determine the "effective" gravitational constant. Thus, we can give a local 
representation of the function (|3*2)l : 

F = 87rG- 1 +C(^ 2 -^)^ ^(1 + 2/(0)) , (35) 

where the dimensionless function y{4>) = 87rG£(</> 2 — 0q) -C 1, at low redshifts (see Thus, on the scales where 

the effective Newtonian constant reduces to the present value of G. (32 ) reduces to (35). 

In this scheme, the present value of the coupling function reduces locally to the inverse of the measured Newtonian 
constant F = (8-kG)^ 1 . 

We will expand equations lj5U|) and (|3*T]) in terms of the perturbative parameter y((po), where the subscript "0" refers 
to the present epoch. Note that 

1 8ttG 



F (1 + y) 



8ttG(1 - y) 



for simplicity, we will assume that the scalar field potential perturbations SV are negligible: in order for this condition 
to hold, even for non-linear scalar field perturbations, the potential must be sufficiently flat. If this is the case, eq. 
(|30|l can be written as 



V 2 $ = 4ttG(1 - y) 



2 A V% 
a dpi 



8ttG 



(36) 



First of all, let us analyze the effect of the background scalar field. In the case of an unperturbed non-minimally 
coupled scalar field, the Poisson equation is modified by a factor (1 — y) with respect to the standard case; this 
correction is scale- independent and proportional to the coupling parameter £ at any redshift, reducing to zero at the 
present time. However, at redshifts relevant for the onset of structure formation, we expect the gravitational force 
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properties to be affected by this correction, which can be positive or negative depending on the sign of the coupling 
parameter £. 

In order to appreciate any substantial correction to the standard Poisson equation which may be induced by the 
perturbations in the coupled scalar field, we have to evaluate under which conditions, and on which length scales L, 
if any, the two terms in square brackets in eq. (|36[1 are comparable. Note that 

AnGa 2 8p m — j{ 2 ^£l!l = L H ~ 2 6 m , 

Pm 

Lh being the Hubble length and 6 m the Dark Matter density contrast. Now, on the scales of clusters, i.e. length-scales 
~ 1CP 3 L#, where matter fluctuations are becoming non-linear today, in order for the two terms on the rhs of eq. 1361) 
to be comparable, it would be sufficient to produce fluctuations Sy ~ 10~ 6 . On galactic scales, overdensities are of 
order 10 4 , the typical scale being 10~ 5 Lh, so again a perturbation Sy ~ 10 -6 would affect at a considerable level the 
Poisson equation. 

However, we are not free to establish the amount of fluctuations in the scalar field component: there are observational 
constraints, coming from the upper limits on the time variation of the gravitational constant and from solar-system 
limits on the JBD parameter ujjbd 0]- The latter is the major constraint, since the Newtonian limit of the Klein- 
Gordon equation is not affected by time variations of the scalar field; we must have 

— !— = ^r = 4C 2 ^8ttG < 2 • 1(T 4 (37) 
where we selected a conservative lower limit lojbd > 5000 = LOi ower . The constraint on ujjbd is thus translated into 



an upper limit on the combination £<^o, since y^jg D = ±2(^00 V87rG; for positive £, it must be 



which gives 



1 / ; — , n \ w, 1 

*o/Wa = g \J^GujBD = C0o8ttG < V8^G 7 • 10~ 3 = 1 (39) 

It is important to emphasize that the condition above has to be thought as applying to the total (background plus 
perturbations) quantities F,F',<p at the present time; i.e. it does not represent a constraint on the present value 
of F^/Fb only, since this value could not be observationally discerned from possible local fluctuations; rather, the 
constraint (|39|l applies to the total ratio Fg/2Fo, which appears in eq. I|31|l . As we will see, this reflects into a 
constraint on £<^o rather than on C4>b,o- 

Hereafter, we will omit the subscript "0" , referring to the value at present (or at low-redshifts) of the quantities 
entering into the Klein-Gordon equation. 

To check whether solar-system constraints on F' /F are compatible with the requirement Sy ~ 10 -6 , thus allowing for 
potentially observable imprints on the structure formation process, we found it convenient to separate the regime of 
linear S<f> ■ V8ttG -C 1 from the more general case 8<f> ■ \/8ttG £ 1. Let us define the dimensionless scalar field variable 
4> = \/8ttG(/). In the first case, inserting the upper limit (|39Jl into eq. (|31|l . we have 



\\f^Lr(^n 2 + 6Cw 2 54>) , (40) 



2 

where the linearity of S<f> allowed to neglect SF with respect to F^ in the factor in front of the brackets, as well as the 
term proportional to \5<p\ 2 ; the last term on the RHS of eq. I|31|) can be neglected on scales much smaller than the 
horizon. The quantities in eq. 140(1 have to be understood as their present values. It is evident that 

s^l\f^L r s m (^-) 2 ; (41) 



L 



H 
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as we have seen, one has a typical value of S m (L/Lh) 2 ~ I0~ 6 on galaxy and cluster scales: therefore, the linear 
value of <5</>v8ttG on those scales is forced to be of the order of 10 -8 (corresponding to Sy ~ 10~ n ), which is orders 
of magnitude below the value required for scalar field perturbations to affect the Poisson equation in a considerable 
way, so that the only effect on the Poisson equation would be given by the overall factor (1 — y), which is active at 
low redshifts even in the absence of field perturbations. 

To analyze the most general case, we will reject the linearity assumption for S<p. On scales much smaller than the 
horizon, and neglecting the field potential, the Klein-Gordon equation 131|) for the model (|35|l reduces in Fourier space 
to 



J JBD 



L x 2 



— ) +3C<ty(2<£- 



(42) 



where we used the non-linear expression for 5F = F — Fb = £8<fi(S(j) + 20b) = (Scj)(2(j> — S(f>). 

Without loss of generality, we can assume that today 4> — 1, because we have only a constraint on the product C,4>; 
with this choice, £ < w ;OT 1 / f 2 „/2 ~ 10~ 3 and eq. l|l2*)l becomes, dividing by \JujJg D ^ 0, 

^JJ^ 54, + 64> 2 + 6 m (J^j ~0 (43) 

On galaxy and clusters scales, a typical value is S m (i^j ~ 10~ 6 , much smaller than the lower limit on 
y/ojjBD ^ \f^iower ~ 10 2 ; the two solutions of l|43|l are therefore 8<p ~ (which trivially corresponds to the so- 
lution previously found under the assumption of linearity), and Sip ~ —^/ojjbd ~ 10 2 , which is a strongly non-linear 
scalar field perturbation; as for the Poisson equation l|36|l. the term Sy is, for S<j> ~ 10 2 , of order 10 , much bigger 
than the value required for the two terms in brackets in eq. i|36|) to be comparable. The fact that such a value does 
not depend on the matter overdensity at any redshift, makes us argue that it is unphysical, and that the Extended 
Quintessence perturbations can only be linear as long as the constraint (|37|l applies. 

However, we want to stress that an initially small scalar field perturbation could, in this model, grow and become 
non-linear on sub-horizon scales, since the limit (|37|l is only restricted to the Solar System neighborhood: no constraints 
on spatial fluctuations of the field are available on larger scales. 

Note again that this kind of "growing" solution is only allowed in the non-minimally coupled case, since we 



assumed y Uj BD 7^ 0: the "Extended Quintessence" model could in principle admit substantial perturbations of the 
scalar field, as their time evolution, according to eq. (|23|l . is sourced by a non- vanishing term which is ultimately 
related to the non- vanishing anisotropic stress (see Appendix 151 and K^i and to the matter perturbations themselves. 



V. CONCLUSIONS 

In this paper we extended the Newtonian approximation to the class of scalar-tensor theories of gravity, where a 
non-minimally coupled scalar field is assumed to be responsible for the cosmic acceleration today. We obtained the 
equations relevant for simulations of gravitational clustering in these cosmological scenarios. The Newtonian Poisson 
equation acquires new contributions from scalar field perturbations, which may turn out to affect the gravitational 
collapse in an interesting way. 

As already argued in a previous study of linear perturbations in these models , a substantial part of scalar field 
perturbations, along with perturbations in its stress-energy tensor components, is powered by perturbations in the 
matter component: this "gravitational dragging" is only due to the coupling of the scalar field with the Ricci scalar 
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(and, ultimately, with the perturbed matter), and can drive scalar field perturbations into the non-linear regime. 
On the contrary, self-interaction originated perturbations in the scalar field will unavoidably damp out on the scales 
relevant for structure formation, as in the case of minimally-coupled fields. 

In order to get specific insights on the role of gravitational dragging, we specialized the final equations to the "Extended 
Quintessence" model of |21| . where the non-minimal coupling strength is quantified by a dimensionless coupling 
parameter. Without adding a potential for the scalar field, we analyzed the result of the Poisson equation for the 
gravitational potentials and for the scalar field perturbations. We isolated and quantified the effect of the time- 
variation of the effective gravitational constant in the equation for clustering produced by the time evolution 
of the background value of the scalar field. We also analyzed the effect of perturbations of the effective gravitational 
constant, finding that, in Extended Quintessence, the scalar field perturbations are prevented from growing non-linear 
by the Solar System constraints on the coupling parameter. Namely, structure formation is locally not affected by the 
modified gravity at a detectable level; in principle, however, scalar-field non linearities are not precluded even in the 
Extended Quintessence scenario, at least on those scales where the coupling parameter is allowed to vary over a less 
restrictive range of values. Further analysis is required in order to quantify the effects on forming structures in more 
general scalar-tensor theories. 
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APPENDIX A: PERTURBATIONS IN ENERGY DENSITY AND PRESSURE 



Let us decompose the scalar field stress-energy tensor as follows |29|: 



T', ■=<>■'<>., --%(<!>,,<'' +2V) : T; ■ : F-[ : T 



iM 



3 



ri (grav) ( 1 . 



(Al) 



The background and perturbed energy densities and isotropic pressure, defined through Tq = — p and T\ — p/3, turn 
out to be 

• Minimal coupling 



(A2) 



Sp<T C) = - 2$ ) + *M) + + 2*)|V<50| 2 + SV - ^ 



2a 2 



(A3) 



S P ( r C) = ^2 (! - 2*)(<ty" + 2fo6</>) - (1 + 2*)|V^| 2 -SV- $^ 



2a 



6a 2 



(A4) 



• Non-minimal coupling 



(nmc) >£ j% (nmc) ^ , tj 



(A5) 



Sp (nmc) = [ 3j p b ^ + 2n <f>) + 3SF{i> + 2H$ -H)- VWSF + (1 + 2^)V 2 5F] (A6) 
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2 1 F 1 2 

Sp (nmc) = — ^ + ^ _ 2$ ^^ _ h^ H _ m ^ _ | _ 2 ^ + V(5FV ( 2 ^ _ 3$ ) _ ( X + 2 *)V 2 ^ (A7) 



Gravitational 



^( r -) =3 ^ (K -i_ jFb) . p^)=-^( /t -i-F b )(2W + W a ) ; (A8) 



^ ( 9 r«,) = _ 3 ^ SF + 1 ( K -i _ Ffa _ ( 5 j f)(2V 2 * - 6H 2 $ - 6ft*) (A9) 



^(f™'") = (2W + ft 2 )^ + - F b - SF) 



(2H + H 2 )<P + + * + + -V 2 ($ - *) 



(A10) 



In the equations above, the subscript 6 refers to background quantities. 
In the Newtonian limit, we have 

(5p^ + 3S P ^ mc ^ -26V (All) 
(dp* + 3<5p )« — (A12) 

(Sp^ + iSp^^^iK-'-F)-^- . (A13) 

At the same time, it is easy to verify that the Newtonian limit of the perturbed Einstein equation 5R® — k5Sq gives 

a 2 

V 2 $ = —K(6p + 3dp) (A14) 

where 8p and Sp include contributions from matter and from scalar field (minimal coupling, non-minimal coupling 
and gravitational parts). Substituting the limits <|A11|) - I|A13(1 into (|A14(i . and comparing with eq. (|30|l we see that 
the contribution —2a 2 8V in eq. (|30|l originates from the minimal coupling, the term — V 2 <5-F originates from the 
non-minimal coupling, while the gravitational contribution can be absorbed in the quantity F . Looking at eq. 
(I3UII . we see that the role of the non-minimal coupling is twofold: on the one hand, it generates the energy and 
pressure contribution — V 2 SF, which is non-vanishing only if the scalar field fluctuations are non vanishing; on the 
other hand, it also generates a modification of the gravitational constant, where the effective gravitational constant 
F^ 1 is generally different from k even if the scalar field does not have fluctuations. 

APPENDIX B: ANISOTROPIC STRESS 

The unperturbed space-space component of the Einstein tensor has only trace component: 

G k--^+l^) S ) (Bl) 
On the other hand, the perturbations SG) have a trace and a traceless part: 



3 [trace) q2 ^3 



(2H + H 2 )<f> + + * + 2Tii> + i(V 2 $ - V 2 *) 



(B2) 
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and 



^^3 (traceless) ^2 



_ 7 ^($_^) jfc + I V 2 ($-*)^ 



(vect. and tensor pert.) (B3) 



Here, we allowed for tensor and vector perturbations, which we do not write explicitly (for details, see 62] ). 
The total stress energy tensor for the scalar field in the Extended Quintessence model can be decomposed as 

rpi_i rpfi(mc) , rpi_i{nmc) . rpfiigrav) 

v v ' v * v 

The anisotropic stress is defined in terms of the spatial components of the stress-energy tensor; the latter can always 
be decomposed as the sum of a trace tensor and a traceless tensor: 

T;= Ph \S}+ir L 5) + ir T )}, (B4) 

where p^ is the fluid pressure in the unperturbed state; itl is the isotropic pressure perturbation (a commonly used 
definition is ttl = Sp/pb), and pb^r) is the anisotropic stress tensor, defined as the traceless part of Tj: 

1 

3' 



We are working in configuration space, so ttl and ttt] are functions of (x, rj). 
Note that in linear theory we can write 

Trrj = {l lk d t d k - |v 2 )^ T 
ttt being the amplitude of the anisotropic stress perturbation, see eg. 0], 

In the following, we will consider linear perturbations in the metric (described in the conformal Newtonian gauge), 
but the scalar field perturbations will generally be non linear. 
For the minimally coupled field, with u> = 1, the tensor Ty- mc ' is given by 

T i = <f> ~ 2 Sj (<t>,c<t>' + 2V) = Tj (trace) + T 3 {traceless) 



where 



and 



ntL) = 3^(1 + 2*)|V^| 2 ^. - \S){^ + 2V) 



Titleless) = j^ 1 + m^d^H - ^(1 + 2n^H\H) 



Note that, since <pb does not depend on the spatial coordinates, di4> = diScj)', for this reason, in the linear regime (i.e., 
for linear perturbations of the scalar field), the anisotropic stress from a minimally coupled scalar field is negligible, 
being second order in the scalar field perturbations. In general, the anisotropic stress from a minimally coupled scalar 
field is of order 0(\VS(j)\) 2 . 

For the non-minimal coupling tensor, 



rpM mC ) _ p,i _ glp,C _ rpi( nmC ) , rpi( nmC ) 

3 '-3 3 - c 3 [trace) 3 (traceless) 



with 



Anrnc) _ 1 j 
3 (trace) ~ 3fl 2 °3 



V^V* - (1 + 2*)V 2 <5F + F(m ~ 3* - 6H$) 
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and, after some algebra, it can be shown that 



T. 



(nmc) 



i 



J (traceless) 



(1 + 2*) 



a- 

(1 + 2*) 



-S^ 2 F 



7^*^+7^**7^ 



i A i 



^ -o'JV.FV* 

3a z J 



1 



y lK yj5F tk - -<S!WFV* 

o J 



We see from here that the anisotropic stress from the non-minimally coupled field contains terms which are proportional 
to the spatial gradients of F; in general, these terms are of the order 0(|V<50|), so the anisotropic stress for a non- 
minimally coupled scalar field would survive even in the case of linear perturbations of the field. 
Finally, the gravitational term is 



T 



(grav) 



1 



- F G\ = T. 



i(grav) 



T 



(grav) 



3 (trace) 3 (traceless) 



Using equations IjBljl . l|B2jl and (|B3|) - we have 



T 



(grav) 



and 



3 (trace) 



i(grav) 



— S 1 

,2 3 



-(H + -H 2 ) + {2H + ft 2 )* + W* + * + 2H* 



r (V 2 *- V 2 *) 



3 {traceless) 





i 


(H 






2 

a 



1 



3 k 



V 2 (*- *)5* 



(vect. and tensor pert.) 



Therefore, the total anisotropic stress, without taking into account vector and tensor perturbations of the metric, 
will be given by the sum of the three contributions from the minimal coupling, the non-minimal coupling and the 
gravitational terms. Since its background value is zero, we have 

1 



ST ^c +nmC+ grav = 1 (] 2$) f ^ ^ , 
3 traceless (2 



3a- 



■(1 + 2*)|V<50| 2 <5?- 



(1 + 2*) 



F 



l lk 5F. kj - ^5 l ^J 2 5F 



1 

3^ 



1 



f^ tj SF k - -(5!V<5FV* 



y fc ($ _ ^) 



3" J 

{vect. and tensor pert.) 



Using eq. I|B3|> , and writing the gravitational contribution on the left-hand side of the traceless part of the perturbed 
space-space Einstein equation 



we obtain: 



3 (traceless) ^ ^^3 (traceless) 



-y fc ($-*) Jfe + -v 2 (*-*) ( 5;- 



= (1 + 2*)7 lfc <9 fc o v , lv 



— (1 + 2*)|V^| 2 ^.+ 



(B5) 



(1 + 2*) 



f k 5F k] - -S)V 2 SF 



j ik ^ d SF k - -jjVtfFV* 



where contributions from vectors and tensors have been omitted. 

As an important consequence, we can note that, in the case of a minimally coupled scalar field (i.e., F = k^ 1 ), one 
has 



* = * 



if 5<fi are linear perturbations in the minimally coupled scalar field, the difference between the two gravitational 
potentials is negligible in a first-order theory. 
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APPENDIX C: NEWTONIAN APPROXIMATION FOR THE MINIMALLY COUPLED SCALAR FIELD 

IN LINEAR PERTURBATION THEORY 

In order to better understand the important role of the non-minimal coupling, let us see what to expect to be the 
perturbation behavior in the case of a minimally coupled scalar field. 

First of all, in the minimally coupled case, the anisotropic pressure perturbation is of order 0(5<fi) 2 , constraining the 
two gravitational potentials $ and if? to differ by terms of the same order. Thus, we expect that the (Newtonian) 
Poisson equations H29fl . (|30|l will be equivalent in linear theory, for a minimally coupled scalar field. 
Taking the trace of the perturbed space-space Einstein equation 5G\ = k6T\ we obtain, for the case of minimally- 
coupled scalar field (k = 8ttGn), 



{2H + H 2 )§ + H§ + "f? + 21ii? + -{V 2 <f> - V 2 *) = - 

3 2 



1(1 - 2$)(<5</) 2 + 20b<ty) - $>4>l - i|V^| 2 - a 2 6V 



(CI) 



Here V 2 (<& — VP) ~ O|V50| 2 can be neglected for linear perturbation of the field. As a consequence, a 2 5V will be 
of order Tt 2 (in eq. IjClfl . 6(f) denotes non-relativistic perturbations); by differentiating a 2 SV with respect to the 
conformal time, one can show that a 2 5V ~ 0(TL 2 ): inserting this result in the Newtonian limit of the perturbed 
Klein-Gordon equation for the non-relativistic scalar field perturbations \2'6\i . we can see that, for linear perturbations 



\7 2 5<f) ~ a 2 5V' ~ 0{H 2 ) , 

which is negligible in the Newtonian limit of scales smaller than the horizon ; thus, the Klein-Gordon equations 
implies that the (minimally coupled) scalar field perturbations, in linear theory, will be negligible on those scales. For 
this reason, the scalar field will behave as a homogeneous component on the scales relevant for structure formation. 
Correspondingly, and the two equations (|50)l and Q29[l. with F — (87rG*) _1 , will be identical (up to terms of 

order {H 2 ,6<p 2 )). 
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